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Vertical  force  acting  on  body 
Arbitrary  function 
Unit  vector  in  x-direction 
Body  length 

Pitch  moment  acting  on  body 

Unit  normal  vector,  positive  into  the  fluid 

Polar  radius  coordinate 
Body  radius  as  function  of  x 


Maximum  body  radius 
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Sectional  area,  S(x)  =  Tr[rQ(x)] 
Forward  velocity  of  the  body 
Velocity  vector  of  the  fluid 
Cartesian  coordinates  (cf.  fig.  1) 
Height  of  body  axis  above  the  wall 
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2D 


Dummy  variable  corresponding  to  the  coordinate 

Polar  angular  coordinate 

Fluid  density 

Velocity  potential 
Two-dimensional  velocity  potential 


ABSTRACT 


Analytical  results  are  derived  for  the  vertical  force  and 
pitch  moment  acting  on  a  body  of  revolution  which  is  moving 
with  constant  forward  velocity  parallel  to  an  infinite  hori¬ 
zontal  wall.  It  is  assumed  that  the  fluid  is  ideal  and  incom¬ 
pressible  and  that  the  body  is  slender  and  situated  close  to 
the  wall.  The  force  and  moment  are  given  in  terms  of  simple 
integrals  involving  the  body  shape  and  separation  distance 
from  the  wall.  Calculations  are  made  for  a  spheroid,  and  it 
is  shown  that  the  vertical  force  rises  rapidly  as  the  body  is 
brought  progressively  closer  to  the  wall. 

INTRODUCTION 


It  is  well  known  in  naval  architecture  that  when  a  body  moves 
through  a  fluid  parallel  to  a  wall,  an  interaction  force  and  moment  will 
act  on  the  body  because  of  the  presence  of  the  wall.  Normally  the  force 
will  be  positive  in  the  direction  toward  the  wall,  due  to  the  venturi  effect 
of  the  accelerated  flow  between  the  body  and  the  wall,  and  this  force  will 
become  more  pronounced  as  the  body  is  moved  closer  to  the  wall.  The  moment 
will  general  1/  increase  also,  under  these  circumstances,  but  its  sign  is 
not  definite  and  will  depend  on  the  detailed  shape  of  the  body;  for  a  body 
with  fore-and-aft  symmetry  the  moment  will  vanish,  under  the  assumptions 
of  ideal  flow,  because  of  the  corresponding  symmetry  of  the  resulting  flow. 

There  are  several  instances  in  which  this  problem  may  be  of 
practical  interest,  particularly  for  a  ship  moving  in  restricted  water  near 
a  bank  or  wall  of  a  canal  (and  where  at  low  speeds  the  effects  of  the  free 
surface  are  not  important) .  Recent  interest  in  oceanographic  submarines 
has  posed  a  similar  problem  for  a  vehicle  moving  in  close  proximity  to  the 
bottom. 

To  obtain  a  complete  analytical  prediction  for  the  force  and  moment 
acting  on  a  body  near  a  wall,  it  is  necessary  to  idealize  the  problem.  With 
this  goal  in  mind  we  shall  assume  not  only  that  the  fluid  is  ideal,  incom¬ 
pressible,  and  unbounded  except  by  the  wall  and  the  body,  but  also  that  the 
body  is  axisymmetric ,  slender,  and  in  close  proximity  to  the  wall.  The 
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last  assumptions  permit  us  to  use  the  s lender-body  approach  of  aero¬ 
dynamics^  to  obtain  an  explicit  solution  of  the  boundary-value  problem  for 
the  velocity  field,  in  terms  of  a  source  distribution  inside  the  body  plus 
a  corresponding  image  source  system  belew  the  wall.  The  strength  of  this 
source  distribution  is  unchanged  by  the  presence  of  the  wall,  but  the 
effect  of  the  wall  is  not  only  to  introduce  the  image  system,  but  also  to 
require  that  the  original  source  distr.bu.ion  be  offset  from  the  body  axis 
toward  the  wall:  this  can  be  inferred  from  the  corresponding  two-dimen¬ 
sional  problem  in  bipolar  coordinates.  If  the  body  were  not  axisymmet r i c , 
a  more  complicated  singularity  distribution  would  be  required,  and  this 
could  only  be  obtained  from  a  numerical  scheme.  Nevertheless,  it  is  felt 
that  the  present  results  should  be  valid  for  bodies  which  do  not  depart 
too  much  from  axial  symmetry,  as  in  the  case  of  a  submarine. 

Having  found  the  velocity  potential,  in  terms  of  a  three-dimen¬ 
sional  source  distribution  within  the  body  and  its  image,  we  can  easily 
find  the  force  and  moment  acting  on  the  body  by  the  use  of  Lagally's 
theorem,2  and  the  final  results  take  the  firm  of  integrals  over  the  body 
length  involving  the  cross-sectional  area  of  the  body  and  its  first  deriv¬ 
ative,  and  the  distance  of  the  body  from  the  wall.  These  results  are 
illustrated  for  a  slender  spheroid. 

Several  years  ago  Eisenberg^  considered  the  same  problem  for  a 
spheroid,  and  although  his  results  were  limited  to  finding  the  pressure 
distribution  on  the  body,  the  force  and  moment  could  be  found  therefrom  by 
integration.  However,  there  is  a  significantly  different  assumption  in 
Eisenberg's  analysis,  from  those  made  here,  .n  that  the  flow  is  approximated 
by  that  due  to  a  single  spheroid  in  an  infinite  fluid  plus  the  corresponding 
disturbance  due  to  an  image  spheroid,  without  regard  for  the  fact  that 
there  will  be  an  induced  velocity  on  each  spheroid  due  to  its  image  which 
must  be  accounted  for.  In  fact  Eisenberg's  approximation  will  be  valid 
provided  the  two  spheroids  are  far  apart,  whereas  ours,  based  upon  the 
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lender-body  theory,  is  valid  provided  the  two  bodies  are  close  together. 

:  is  felt  that  the  latter  is  more  appropriate  in  the  present  context, 
nee  the  force  and  moment  under  consideration  increase  rapidly  as  the 
>dy  is  brought  closer  to  the  wall  and  are  comparatively  negligible  if 
te  body  axis  is  greater  than  one  or  two  diameters  from  the  wall. 

THE  VELOCITY  POTENTIAL 

We  shall  employ  Cartesian  coordinates  (x,y,z),  with  z  =  0  the  plane 
the  wall  and  with  the  body  axis  on  tht.  line  y  =  0,  z  =  Zq.  Thus  the 
dv  axis  is  parallel  to  the  x-axis  and  is  a  distance  z^  from  the  wall 
igure  1) .  The  body  is  assumed  to  be  defined  by  a  smooth  continuous 
dius  function  rQ(x) ,  according  to  the  equation 

rQ(x)  *  [y  +  (z  -  zo)  ]  ,  for  - y <  x  <  y  ,  [1] 


Figure  1  -  The  Coordinate  System  and  Flow  Geometry 

re  L  is  the  body  length  and  by  assumption  rQ(±  =  o.  Because  of  the 
nder  body  assumption  and  because  the  wall  is  in  the  near  field, 

L  <<  1  and  z  /L  <<  1. 


The  body  is  considered  to  be  fixed  in  space,  with  a  free  stream 
of  velocity  U  flowing  past  it  in  the  -x  direction.  The  velocity  field  is 
defined  by  the  vector 


V  =  V  *(x,y,z)  -  Ui  [2] 

where  <p  is  the  velocity  potential  of  the  perturbation  flow.  This  potential 
must  satisfy  Laplace's  equation,  the  boundary  condition 

=0  on  the  wall  :  *  0, 

a  Z 

and  the  boundary  condition 

U  cos(n,x)  =  -  U  r  '  (x) 

3n  o 


on  the  body.  Here  n  is  the  unit  normal  vector  out  of  the  body  and  r  ' (x)  3 

1  o 

drQ/dx.  In  addition,  must  vanish  at  infinity.  In  the  near  field,  close 
to  the  body,  it  follows  from  slender  body  theory*  that 


<t>(x,y,z)  =  <t>2D(y,z;x)  +  f(x) ,  [3] 

where  is  the  potential  of  a  two-dimensional  flow  in  the  "cross  -  flow" 
plane.  The  function  f(x)  can  be  found,  but  since  it  depends  only  on  x  it 
does  not  contribute  to  the  vertical  force  or  to  the  pitch  moment. 

The  potential  must  satisfy  the  body  boundary  condition 


3  $ 


2D 


3n 


-  U  r  '  r x) 
o 


[J] 


Physically  this  is  the  velocity  potential  for  the  flow  of  fluid  out  (or  in! 

2  ?  2 

through  the  periphery  of  the  circle  y  +  (z  -  =  rQ  ,  in  the  presence 

of  a  wall  at  z  =  0.  By  reflection  this  is  equivalent  to  the  flow  out  of 

(or  in  through)  two  circles  with  centers  at  z  =  -z  .  To  solve  this  two- 

o  * 

dimensional  problem,  we  look  for  a  solution  of  the  form 


Some  motivation  for  attempting  to  find  a  solution  of  this  form  can  be 
fc  md  from  the  properties  of  bipolar  coordinate  systems. 
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4>2q  (y,z;x)  =  C  log  |  [y2  +  (z-a)2]  [y2  +  (z  ♦  a)2]  ?  [5] 

lere  C  and  a  are  constants  (which  may  depend  on  x)  .  Physically  this 
)tential  represents  the  flow  due  to  a  source  at  y  =  0,  2  =  a,  plus  an 
nage  source  at  y  =  0,  2  =  -  a.  Changing  to  polar  coordinates  (r,  9),  with 
le  origin  at  the  center  of  tfie  original  circle,  gives 

y  =  r  sin  9 

2  =  2  +  r  cos  9 , 

o  ’ 


1  it  follows  that 


2  2  2  2 
4>2p  ( r ,  9 ;  x)  =  C  log  (r  +  2  z  r  cos  9  +  zq  *  a  ) 


1  2  ,  .2 

-  4  a  (2  +  r  cos  9  ) 

o 


can  be  easily  verified  that,  on  the  circle  r  =  rQ(x), 


If  *2D  '  2  C/V 

ovided  the  parameter  a  is  chosen  to  satisfy  the  relation 

2  2  2 
a  =  z  -  r 
o  o 

om  the  boundary  condition  [4]  it  follows  that 


C  =  -  i-  U  r  (x)  r  '  ix)  =  -  f*U  S'(x), 
2  0  0  4  it 


ere  S(x)  =  irro"(x)  is  the  cross-sectional  area  and  S'(x)  =  d  S(x)/dx. 
us  the  desired  two-dimensional  potential  is 

(y,z;x)  =  -  *2“U  S '  (x)  log  {  [y2  ♦  ( 2  -  a) 2  ]  [y2  *(z  *  a)2]]  [7] 


a(x)  =  z2  -  [rQ(x) ] 2 
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The  three-dimensional  source  distribution  which  is  associated  with 


<t>(x,y,z) 


I  (x-n2 


2 

♦  v  “  +  ( 


[8] 


♦  [ ( x- c J -  *  y“  ♦ 


1  - 1  n 
(z*a)“]  W ‘ 


dr 


where  in  the  integral  the  parameter  a  is  a  function  of  r, .  It  can  be  shown 
that  this  three-dimensional  potential  reduces  to  as  given  by  liquation 
[7],  plus  a  function  only  of  x,  for  points  on  or  near  the  body  surface. 


THE  VERTICAL  EORCE  AND  PITCH  MOMENT 


It  is  now  a  straightforward  matter  to  find  the  vertical  force  F, 
and  pitch  moment  M  from  Lagally's  theorem.”  For  this  purpose  we  note 
from  [8]  that  the  body  is  generated  by  a  distribution  of  three-dimensional 
sources,  of  strength 

7;  U  S'(x), 

in  the  presence  of  an  "external  flow"  due  to  the  free  stream  of  velocity 
-  Ui ,  where  i  is  a  unit  vector  in  the  x-direction,  and  an  additional  ex¬ 
ternal  flow  due  to  the  image  potential 


S'(£)  [(x-r,)”  +  y"  +  (z  +  a)”]  dC- 

~2~ 

From  Lagally's  theorem  the  vertical  force  acting  upon  the  body  is  the 
integral  of  the  differential  force  on  each  source  element,  which  in  turn 
is  equal  to  4*3  times  the  product  of  the  source  strength  and  the  vctical 
component  of  the  external  flow.  Thus 

a  ,  -1/2  I 

S’t  •)  [(x-fj"  ♦  (:  ♦  a)"]  dU  dx.  [«J 

:  =  a  ( x ) 
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The  quantity  in  brackets  can  be  evaluated  t'or  a  slender  body  since  it  is 
a  vertical  dipole  distribution: 


i-  C2  s. 

32Li 

2 


- 1/2 

S'  (O  [(x-a2  +  (2  +  Ca(C))2]  dc 


2=a(x) 


L 

(O  [a(C)  +  a(x)  ]  [(x-O2  +  (a(£)  *  a(x))2]  d; 


"-J  I3’ 

=  -  2  a(x)  S 


-o/2 


■3/2 

[(x-C)"  ♦  4  (a(x))‘J  d; 


=  -  S'  (x)  /  a(x) . 


Thus,  from  [9],  the  vertical  force  is  given  by  the  integral 

L 


Fz  =  -  4 


47  plj2  f 2  Is’ 

J  L 


(X)  ]■ 


2o2  -  t^*)]2 


-V2 


dx, 


or 


hz  -  ffplj2  J  L  [r0(x)  ro'(x)]2{ 


2  2 
zo  -  [ro(x)J  dx. 


[10] 


Similarly,  for  the  moment, 

L 


M 


->  r  2 

=  noU“^[ro 


2  f  2  2I ”^2 

(x)  rQ'(x)]2  zq2  -  [rQ(x)]M  x  dx. 


[11] 


These  two  integrals  are  the  principal  results  of  the  analysis.  We 
note  that  for  a  body  with  longitudinal  (fore-and-aft)  symmetry  the  integral 
for  the  pitch  moment  M  vanishes,  as  a  consequence  of  the  symmetrical  flow 
vhich  wil  then  occur.  Moreover,  the  vertical  force  r  is  ‘always  less  than 
zero,  in  keeping  with  the  fact  that  the  flow  velocity  will  be  higher 
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between  the  wall  and  the  body,  as  compared  with  the  corrcspc  iding  value 
on  the  "outside"  portion  of  t lie  body.  Thus  there  will  always  be  a  suction 
force  tending  to  attract  the  body  toward  the  wall  and  this  force  will  in¬ 
crease  monotonical ly  to  a  finite  limit  as  :  •»  r  ,  or  the  bodv  is  brought 
in  to  the  wall.  We  can  conclude  from  the  above  integrals  that  for  geo¬ 
metrically  affine  bodies,  having  the  same  cross-sectional  distributions 
but  different  lengths,  the  force  will  be  inversely  proportional  to  the 
length  whereas  the  moment  will  be  independent  of  the  length.  Finally,  we 
note  that  for  bodies  with  flat  ends,  such  as  a  right  circular  cylinder, 
the  above  integrals  are  divergent  as  a  consequence  of  the  "end  effect" 
error  in  the  s lender-body  theory. 

The  pitch  moment  is  more  difficult  to  discuss  qualitatively  since 
the  integrand  of  [11]  is  not  positive-definite,  and  the  sign  of  the 
integral  will  depend  critically  on  the  distribution  of  the  sectional  area 
S(x) ,  and  possibly  also  on  the  separation  distance  :  .  However,  it  seems 
likely  that  for  bodies  with  full  or  blunt  ends  this  integral  will  be 
dominated  by  the  large  values  of  S'(x)  at  the  two  ends,  and  thus  that  the 
pitch  moment  will  be  positive  (bow  down)  if  the  bow  is  blunt  compared  to 
the  stern,  and  vice  versa.  (An  analogous  argument  is  common  in  naval 
architecture  to  explain  the  tendency  of  ships  with  full  sterns  to  tend  to 
head  away  from  a  bank  when  moving  in  restricted  waters.) 

THE  VERTICAL  FORCE  ON  A  SF  1ERI0D 


As  r.n  illustration  we  shall  consider  the  body  to  be  a  slender 

spheroid,  of  maximum  radius  r  ,  so  that 
r  *  om’ 


r0(x) 


r 

om 


(1  -  4x2/L2) 


1/2 


For  convenience  we  introduce  the  nondimensional  parameter 


representing  the  nondimensionalized  vertical  height  of  the  body  axis  above 
the  bottom.  Substituting  in  [10],  we  obtain  the  integral 
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